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Abstract. Percolation theory is usually applied to lattices with a uniform probability 
p that a site is occupied or that a bond is closed. The more general case, where p is a 
function of the position x, has received less attention. Previous studies with long-range 
spatial variations in p{x) have only investigated cases where p has a finite, non-zero 
. gradient at the critical point p c . Here we extend the theory to two-dimensional cases 

in which the gradient can change from zero to infinity. We present scaling laws for the 
width and length of the hull (i.e. the boundary of the spanning cluster). We show that 
the scaling exponents for the width and the length depend on the shape of p(x), but 
they always have a constant ratio 4/3 so that the hull's fractal dimension D = 7/4 is 
qq ■ invariant. On this basis, we derive and verify numerically an asymptotic expression 

| for the probability h(x) that a site at a given distance x from p c is on the hull. 
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Figure 1. Site percolation on a triangular lattice where occupied sites (grey 
hexagons) are more abundant towards increasing x-coordinates. The hull (black curve) 
marks the boundary between the vacant spanning cluster on the left and the occupied 
spanning cluster on the right. The hull width w, defined in ([T]), is the standard 
deviation in the hull's x-coordinates. 



1. Introduction 

Percolation theory is a powerful tool in describing systems of randomly located objects 
that are placed on the sites or bonds of a network [1]. The sites or bonds are 
independently occupied with a probability p or remain vacant with probability 1 — p. 
The task of percolation theory is to identify the properties of the clusters formed by the 
connected components. In most theoretical studies, p is assumed to be equal everywhere 
(i.e. the system is assumed to be uniform). However, in many real- world applications, it 
is more natural to allow some variation in p among the sites or bonds. In spatial systems, 
p often depends on the position. For example, if occupied sites represent particles 
that diffuse from the right-hand edge of a two-dimensional lattice and are absorbed on 
the left-hand edge, then the concentration of particles exhibits a gradient |2J and the 
system looks qualitatively like the lattice in figure [TJ Gradients also frequently arise in 
the spatial distribution of biological populations because the environmental conditions 
may change gradually in space (e.g. from lower to higher elevation along a hillside). If 
occupied sites represent the presence of a species, the population can exhibit a spatial 
transition from isolated individuals to a large connected component 0, 0] ■ 

Given such a transition in connectivity, the focus of gradient percolation is on the 
hull 0|, defined as the interface between adjacent occupied and vacant clusters that both 
fully span the lattice's y-direction (black curve in figure [1]) The hull has been used in 



J If the length L y of the lattice in the y-direction is finite, then there is a finite probability to find more 
than one pair of adjacent occupied and vacant spanning clusters and hence multiple hulls. However, 
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practice, for example, in studies of the surface structure of polymers of the treeline 
dividing woodland from grassland E 8, of the boundary of urban settlements | and 
as an example of a non-trap ping self-avoiding random walk whose properties have been 



>e p 

investigated theoretically |9|, [10[ and with computer simulations 111, [12] . Apart from its 
many applications, gradient percolation is also an active field of theoretical investigation. 
It allows highly accurate measurements of the smallest occupancy probability p c for 
which an infinite spanning cluster of occupied sites exists. The value of p c is equal 
for gradient and non-gradient (i.e. uniform) lattices. Introducing a gradient has the 
advantage that p does not need to be carefully fine-tuned to achieve excellent numerical 



results 13j,ll4|. Furthermore, the geometry of the hull is particularly intriguing. Suppose 
the hull connects the sequence of coordinates (x 1 ,y 1 = 0), . . . , (xi,yi = L y ) delineating 
the boundary of the spanning cluster, where L y is the system size in the ^/-direction. If 
p can be approximated near p c by a linear function p(x) = ax + p c , then the hull width 
w, defined as the standard deviation of the hull's x-coordinates, 



is known to scale as a -4 / 7 in the limit a — V 0. The hull length I is proportional to 
L y a~ 3 / 7 and, on length scales below w, the hull's fractal dimension is D = 7/4 j^]. 



These non-trivial scaling laws - recently proved rigorously by Nolin 15|, |l6| - have 
sparked a lot of interest (see for a survey of the literature) and some researchers 



have generalized the results to gradients in more complicated lattice models [18H23 
So far, all studies have assumed that p(x) can be treated, at least locally, as a linear 
function (figure [2](a)). However, this is frequently not the case in reality. As a -> 0, 
higher order terms in the Taylor expansion of p may become dominant so that w and I 
do not diverge, although this is predicted when interpreting w oc a _4//7 and I oc L y a~ 3 ^ 7 
literally. Local linearization is not possible either if the gradient does not exist at all, 
for instance if p jumps discontinuously from one value to another. 

Here we examine representative cases in which no linear approximation for p(x) 
exists at p c . We investigate the class of functions 

( max(0, p c - a(-x) b ) if x < 0, 

P( x ) = ) ■ , b , . f . n ( 2 ) 

I mm(l, p c + ax ) if x > 

for b > and a > 0. If b = 0, p{x) makes a discontinuous jump over the percolation 
threshold p c (figure [2]^b))J§| For < b < 1, p(x) is continuous, but the gradient at p c 



this probability goes rapidly to zero as L y — > oo [6( , see [1 51 ] for a rigorous proof. We never encountered 
this situation in our numerical simulations with L y > 4096. 

§ For functions defined only on the discrete coordinates of a lattice, we cannot, strictly speaking, 
distinguish continuous from discontinuous behaviour. However, for p{x) defined in ([2]), we can apply 
the following criterion. Let x_ and x+ be the neighbouring left and right coordinates of x in the lattice. 
If and only if p(x) — > p{x-) and p(x) —> p(x+) in the limit a — > + , then the underlying function p 
of real-valued coordinates must be continuous at x. In this case, we also call the discretized function 
continuous. 
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Figure 2. The studied occupancy probabilities p(x). (a): Linear function. This 
has been, at least as a local approximation, the assumption in all previous studies 
of gradient percolation, (b): Discontinuous function (6 = in @). (c) and (d): 
Continuous function with infinite or zero slope at x = (0 < 6 < 1 or i > 1, 
respectively) . 



is infinite (figure EJc)). If 6 > 1, the gradient is zero (figure Efd)), so that \p(x) — p c \ 
approaches zero faster than any linear function at x = 0. We have chosen the coefficient 
a in to be equal on the positive and negative x-axis. This is a convenient choice 
because, as we show below, the probability to find the hull at a coordinate x is in this 
case symmetric about x — 0. The results described below can be easily generalized to 
cases where the coefficients differ in the positive and negative direction. We present 
numerical simulations for site percolation on triangular lattices, where p c = 1/2, and 
apply periodic boundary conditions in the ^-direction (L y > 4096). For a comparison, 
we also performed simulations for site and bond percolation on square lattices. We 
found the same scaling exponents, indicating that they are independent of the lattice 



type and presumably also valid in continuum models 24 



2. Hull width, length and fractal dimension 

Let us first consider the hull width w in the discontinuous case (6 = in (j2J)). Because 
a > 0, the hull must navigate around those occupied clusters to the left that are not part 
of the occupied spanning cluster (figure [3]). Similarly, on the right-hand side, excursions 
of the hull must steer clear of smaller vacant islands that have not merged with the 
vacant spanning cluster. The typical linear size of these smaller clusters is given by the 
correlation length £ which scales in uniform systems oc \p — Vr\~ v if p ~ p c - Because 
the critical exponent v is equal to 4/3 in two dimensions [l|, 0] and because the hull 
width scales in proportion to the cluster sizes, we expect w oc a~ 4 / 3 . Our numerical 
simulations confirm this intuition (black filled circles in figure Hk). 

If b > 0, the situation is more complicated because p(x), and hence also £(p(x)), 
change gradually as a function of x. For 6 = 1, Sapoval et al. [2] have pointed out that 
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Figure 3. Illustration of the relation between hull width w and correlation length 
£ near a; = 0, where the hull meanders between islands of vacant and occupied sites. 
The linear size of these obstacles is of the order of £. Therefore, the hull width w, 
which measures the fluctuations of the hull around x = 0, is also approximately £. 



we can still derive the correct scaling behaviour by assuming that w is proportional 
to the correlation length at w. The self-consistent relation w oc £(p(w)), where the 
local correlation length satisfies £(p(x)) oc \p(x) — p c \ ~ v and p{x) = p c + ax, yields the 
aforementioned w oc a~ v /^ = a" 4 / 7 (black open squares in figure Hh). Carrying out 
the calculation more generally for p(x) given by ([2]) with arbitrary b and assuming a to 
be small, we obtain 

w(xa-^ ub+1 \ (3) 

The same self-consistency requirement has been noticed and tested for the Ising quantum 



chain [28|, |29[ , where the exponent v however has a value different from the exponent 
in percolation. In figure H^a), we compare data from numerical simulations with the 
prediction of equation ([3]) based on the two-dimensional percolation exponent v = 4/3. 
The data show excellent agreement. 

There is a similar scaling relation between the normalized hull length l/L y and 
a, where I is the number of steps in the hull in a lattice of vertical dimension L y . 
The scaling of l/L y cannot be derived as intuitively as the scaling of the hull width. 
However, previous numerical evidence in the case 6=1 has supported the hypothesis 



l/Ly oc a 3 / 7 |2j, now known to be correct 15|. For arbitrary b, our data (figure Hb) 
indicate the more general relation 

l/L y oca^ b+1 \ (4) 

which includes the scaling law for b = 1 special case. 

The scaling of I and w is related to the fractal dimension D of the hull as follows. 
For a given value a = a 1 and system size L y> i, let us denote the hull length by li and the 
width by Wi. If we coarse-grain the hull by measuring it with rulers whose length is A 
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Figure 4. (a) The hull width w versus a for different values of b in ^ from computer 
simulations (shown as circles, squares, and diamonds). The lines indicate the expected 
power law w oc a~ B with B — vj(yb + 1). The number next to each line states the 
value of B. Error bars are smaller than the symbol sizes, (b) The normalized hull 
length l/Ly scales oc a~ c with C = l/(i/b+ 1) as displayed by the lines, (c) The 
fractal dimension D of the hull is determined with the equipaced polygon method 25 1 
described in the appendix. In this method, we compute the mean Euclidean distance 
d(k) between hull coordinates that are k steps apart. A numerical estimate of D can 
be obtained by fitting d(k) oc k 1 ^ . We observe that D — 7/4 is independent of b. 
Inset: D — 1.75, numerically calculated as the derivative d(logfc)/d(logd), is within 
the error bars for almost three orders of magnitude of k. 
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times longer than the lattice constant, we need I2 oc A D li rulers to cover the full length 



of the hull 30]. The width and the system size, however, are linear objects: it requires 
w 2 oc A^wi rulers to span across the hull's x- and L y ^ 2 oc A~ x L y ^ rulers to cover its 
^/-extent. Thus, measured in units that are increased by a factor A, the hull length, 
width and system size all appear smaller by either a factor A~ D or A~ x . Now suppose 
that the hull geometry scales as w oc a~ B and I oc L y a~ c . We obtain the same w 2 and 
Ly j2 as before if the lattice constant remains our fundamental length scale, but we first 
replace a\ by a 2 oc A x l B a\ and then, from the original system of length L y> \, we cut 
out a strip of size A~ l L y ^. In this geometry l 2 oc (A~ 1 L y> x)(A~° /B a^ c ) oc A~ (1+c/jB) /i. 
Comparing this with our previous expression for l 2 , we conclude that 

= ! + ( 5 ) 

Inserting B = vjivb + 1) and C = \j[yb + 1) from equations ([3]) and (J3J) we find 
D = 1 + 1/u = 7/4, independent of b. We also measured the fractal dimension directly 
using an independent numerical method (following [l2, 25] and outlined in the appendix) 
and again find good agreement with D = 7/4 (figureHb). The fractal dimension of cluster 
hulls has recently received a lot of interest. In uniform percolation, where p = p c at all 
sites (i.e. a = in ([2])), D = 7/4 is now known to be exact, thanks to mathematical 
advances in stochastic Loewner evolution [26|. Recent work by Nolin and Werner 27 
proves that the Hausdorff dimension remains 7/4 if p is near, but not exactly equal to 
p c . In constant gradients {b = 1 in (T5])), D = 7/4 was already conjectured in early 
studies I2J and related to the scaling of w and I. Our results demonstrate that D = 7/4 
holds much more generally. Although B = vj[yb + 1) and C = \ j[yb + 1) depend on 
b and thus differ in the different scenarios depicted in figure [2], the ratio B/C = v and 
hence the fractal dimension D = 1 + 1/u are invariant properties of the hull geometry. 



3. Hull density profile 

The hull width and length are simple measures of the hull geometry, but the distribution 
of the hull's x-coordinates contains yet more information. For fixed x, we define the hull 
density h(x) as the fraction of coordinates (2, y) in the lattice that are on the hull. We 
focus first on the case 6 = where p(x) = p c ^f a provided that a is sufficiently small. 
The =F sign is to be interpreted as a minus (plus) sign for negative (positive) x. Our 
simulations for 6 = show that h(x) has a maximum at x = and falls off exponentially 
and symmetrically to both sides (figure [5^) , 

h(x) = EQa u ~ l exp(— Fa u \x\), (6) 

where E and F are positive constants. The exponential decay of h(x) is reminiscent of 
the off-critical correlation function g c {r) oc exp(— |r|/£) in uniform percolation, where 
<7c( r ) denotes the probability that r and the origin are part of the same finite cluster. The 
decay rate {Fay) -1 has the same scaling as the correlation length £, namely oc \p — p c \~ u 
for both p < p c and p > p c . However, there is one important difference between (Fa u ) -1 
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Figure 5. (a) Hull density h as a function of position x for b = in equation 
©. Inset: data collapse for h(x) = a^- v ^ b+ ^h(a- u ^ vb+ ^x). (b) Hull density for 
a piecewise continuous function p{x) with multiple discontinuities. The local value of 
p{x) is given by the grey numbers under the graph. In every interval of constant p, 
h is approximately exponential, resulting in straight lines in this logarithmic plot, (c) 
Hull density for 6 = 4. 
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and the correlation length. In uniform percolation, £ has different amplitudes on the 
sub- and supercritical branch [31] . By contrast, we observe that F is equal for x < 
and x > 0. 

Equation (jSj) has several interesting consequences. First, h(x) is a symmetric 
function (i.e. = h(—x)). We observed this symmetry in all investigated lattice 

types (triangular site percolation, square site and bond percolation). Second, (JBJ) is 
consistent with the scaling relations ([3]) and (j3J) of the width and length for 6 = 0. This 
can be proved by inserting (jHJ) into I = f^° oo h(x)dx and w 2 = (x 2 h(x)dx)/l. Third, 
after appropriate rescaling of h, the function h(x) = a l ~ v h{a~ v x) is independent of a 
(see the data collapse in the inset of figure ^ a)). 

The case b > 0, where p(x) changes continuously, requires a little more effort than 
b = 0, where p(x) is constant except for one discontinuity. Our goal is to approximate 
h(x) for b > by discretizing p(x). To this end, let us consider the intermediate case 
where p has more than one discontinuity: p(x) = for x < x±, p(x) = H\ for x G [xx, X2), 
. . ., p(x) = 7r rt _i for x G [x n -i,x n ) and p(x) = 1 for x n < x with < 7Ti < . . . 7r„_i < 1. 
Let us assume that the length (xj — Xj_i) of each plateau is larger than the exponential 
decay length (F\iTi — p c \ v )~ l - in this case, h consists of piecewise exponential functions 
cx exp(±F|7Tj— p c \ u x) that are continuously joined together at the breakpoints x\, . . . , x n 
(an example is shown in figure [5b). In terms of a differential equation, we can express 
this as 

^-\ogh(x) = ±F\p(x)-p c \\ (7) 
dx 

where F is a constant. The plus sign holds for p(x) < p c ; otherwise the minus sign 
applies. 

We now make two approximations to obtain an analytic expression for h(x) if p(x) 
changes continuously. First we assume that the rate of change in p(x) is small compared 
to the characteristic decay length (F\p(x) — p c \ y )~ l for all x. This allows us to insert 
P directly into ©. We obtain 

^- log h( x) = ±Fa u (TxY b (8) 
dx 

for x in the interval [— (p c /a) l l b , ((1 — p c )/a) 1//fe ]. We know from ()3]) that the hull width 
is a vanishing fraction of the width of this interval for a — > + . Thus, our second 
approximation consists of extending flSJ) over the entire real line. The solution is 

h(x) = E b a {u - l ^^ b+ ^ exp ( — a v \x\ vb+l \ , (9) 

V vb + 1 / 

where we have factored out a^ v ~ 1 ^^ vb+1 ^ from the integration constant so that, according 
to (jl]), the remaining factor E h does not depend on a. 



4. Discussion 



Straightforward integration proves that (Q is consistent with the scaling of the hull 
width in (J3]). Because (EJ) also satisfies the scaling of the length 01]) , it implies the 
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Figure 6. The transformation from h to fb according to (|T0|) shows that the hull 
density decays as exp(— Fa v x wh+1 / {vb + 1)) with the same factor F ps 4.3 for all b. 



correct fractal dimension 7/4 according to fl5]). In addition to the scaling laws, (Q also 
provides an analytic expression for the functional form of the hull density profile: h(x) 



is, in this approximation, proportional to a generalized normal distribution [32[. For 
b = we predict a Gaussian profile. Another special case is b = 0, where we retrieve 
the Laplace distribution of (EJ. For < b < \jv the distribution is leptokurtic (i.e. 
it has a narrower peak and fatter tails than a Gaussian); for b > 1/v it is platykurtic 
(broader peak, thinner tails). 

For any particular b in (Q, the rescaled hull density h(x) = a^ l ~ u ''^ ub+1 'h{a~ v '^ vh+l ^x) 
is independent of a (insets in figures|5]^a) and[5(c)). Rescaling can even go further. Above 
we have assumed that the constant F in (jj]) is independent of b so that the functions 

f b (x) = \og{a {l ~ u)/{vb+l) h[{{ub + l)a- v x) 1 ' { - vh+ % (10) 

ought to be linear functions with the same slope — F for all b. Figure |6] demonstrates 
that this is indeed true for x > 1, which is valid if x > w. From linear regression to f b , 
we estimate F = 4.3 ± 0.1 J] 

For a simple line in figure [6] fits the data very well for all x > and this is 

presumably exact for b = 0. However, as b becomes larger, there are visible departures 
from the asymptotic line for i C 1. In this regime, p(x) varies significantly over length 
scales comparable to the correlation length, thus violating the assumption behind (]8]). 
If 6 ^ 1, we are effectively dealing with a system at the critical point confined to the 
strip —(pc/a) l l b < x < ((1 — f} r ) I a ) l ^ b ■ Some properties of critical cluster geometry in 
a strip are known exactly 33N35|. From the known results for the "one-pinch point 
density" in rectangular domains j36|, it is plausible that h(x) is analytic and has a 
Taylor expansion around x = with leading terms h(x) = h(0) + h"(0)x 2 /2 + O(x A ) and 
h"(0) < 0. As a consequence, as x — > we expect log(h) to be approximately quadratic 



(i.e. log(h(x)) ~ G — H x 2 ), in contrast to log(h(x)) ~ — H 



for \x\ — > oo 



|| F depends on the lattice type, unlike B, C and D. The stated value is for triangular site percolation. 
For square site percolation, we find F = 3.9 ± 0.1, for square bond percolation F = 5.4 ± 0.1. 
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Figure 7. If the hull density is locally approximated as a function proportional to 
exp(— H \x\ J ), the best fitting exponent J increases as vb + 1 in the tail (i.e. as x — > oo), 
but is close to 2 for b > 3/4 in the core (i.e. as x — > 0). 



as predicted by (jHJ). 

We test this hypothesis with nonlinear least-squares regressions of the general form 
G — H\x\ J to fit either the core or the tail of \og(h(x)). The best fitting exponents for 
| a; | — > oo are indeed consistent with J = vb + 1 (figure [7]) for all 6. Performing the 
regression around x = 0, on the other hand, we find a crossover from J = vb + 1 to 
J = 2at6~ l/V = 3/4. This suggests that for b < l/v, (J9]) approximates the hull 
density well, but for b > 1/V we should replace it with 

F 



h(x) E b a 



{u-l)/{vb+l) 



exp 



-a \x 



ub+l 



2v/(vb+l) x 2 



11 



i/6 + l 

This approximation has the observed asymptotic behaviour for both x = and \x\ — > oo: 
h(x) is Gaussian in the center, but drops more rapidly in the tails. In the inset of figure 
[5(c) we have used F = 4.3 from above and fitted Ef, and Kj, to the hull density for 
6 = 4. The regression curve indeed fits the core of the distribution better than (Q and 
also provides an excellent approximation in the tails. 



5. Conclusion 



In summary, we have investigated a generalization of gradient percolation where the 
occupancy probability p(x) is nonlinear at the critical point p c . If, to lowest-order, 
\p(x) — p c \ = a\x\ b for a > and 6 > 0, then the hull width and length scale as 
w oc a~ u ^ ub+1 ^ and I oc a^ l ^ ub+1 \ The hull's fractal dimension D = 7/4 is independent 
of 6. The hull density is symmetric and drops for x far away from p c in proportion to 
exp(— Fa v \x\ ub+l / (vb + 1)), where F is a lattice-dependent constant independent of a 
and 6. Closer to p c , the hull density becomes approximately Gaussian if 6 > 1/V. 
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Appendix. Equipaced polygon method to determine the fractal dimension 

There are many different ways to determine the fractal dimension of the hull. We tested 



several of them and found that the equipaced polygon method of [25| provides a good 
trade-off between simplicity and consistent numerical results. We begin by labeling all 
hull coordinates (x^yi) with a subscript i equal to the number of steps required to 
reach this position along a walk on the hull starting at the bottom of the lattice. At this 
stage, the hull is described by a sequence of coordinates (x\,yi = 0) . . . (xi,yi = L y ), 
where I is the hull length. For a given k and for 1 < i < I — k, we then calculate 
the Euclidean distances di{k) between (xi,yi) and (x i+ k,yi+k)- Finally we compute the 
mean d(k) =_52idi(k) / '(I — k + 1). If the hull is a fractal of dimension D, we expect 
d(k) oc k 1/D yr 

References 

[1] Stauffer D and Aharony A 1991 Introduction to Percolation Theory (London: Taylor 
& Francis) 

[2] Sapoval B, Rosso M and Gouyet J F 1985 J. Physique Lett. 46 L149-56 

[3] Milne B T, Johnson A R, Keitt T H, Hatfield C A, David J and Hraber P T 1996 
Ecology 77 805-821 

[4] Gastner M T, Oborny B, Zimmermann D K and Pruessner G 2009 Am. Nat. 174 
E23-39 

[5] Voss R F 1984 J. Phys. A 17 L373-L377 

[6] Grassberger P 2002 Comput. Phys. Commun. 147 64-70 

[7] Wool R P and Long J M 1993 Macromolecules 26 5227-39 

[8] Makse H A, Andrade J S Jr, Batty M, Havlin S and Stanley H E 1998 Phys. Rev. 
E 58 7054-62 

[9] Weinrib A and Trugman S A 1985 Phys. Rev. B 31 2993-7 
[10] Saleur H and Duplantier B 1987 Phys. Rev. Lett. 58 2325-2328 
[11] Ziff R M, Cummings P T and Stell G 1984 J. Phys. A 17 3009-3017 
[12] Ziff R M 1986 Phys. Rev. Lett. 56 545-548 

[13] Rosso M, Gouyet J F and Sapoval B 1985 Phys. Rev. B 32 6053-4 

[14] Ziff R M and Sapoval B 1986 J. Phys. A 19 L1169-72 

[15] Nolin P 2008 Ann. Probab. 36, 1748-1776 

[16] Nolin P 2008 Electron. J. Probab. 13, 1562-1623 



REFERENCES 13 



[17 

[is; 

[19 
[20 
[21 
[22 
[23 

[24 
[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 
[33 
[34; 
[35 
[36 



Gouyet J-F and Rosso M 2005 Physica A 357 86-96 

Kolb M, Gouyet J F and Sapoval B 1987 Europhys. Lett. 3 33-8 

Boissin N and Herrmann H J 1991 J. Phys. A 24 L43-7 

Lemarchand A, Nainville I and Mareschal M 1996 Europhys. Lett. 36 227-31 

Sapoval B, Baldassarri A and Gabrielli A 2004 Phys. Rev. Lett. 93 098501 

Loscar E S, Guisoni N and Albano E V 2009 Phys. Rev. E 80 051123 

Gastner M T, Oborny B, Ryabov A B and Blasius B 2011 Phys. Rev. Lett. 106 

128103 

Quintanilla J, Torquato S and Ziff R M 2000 J. Phys. A 33 L399-L407 

Kaye B H and Clark G G 1985 Part. Part. Syst. Char. 2 143-8 

Smirnov S and Werner W 2001 Math. Res. Lett. 8 729-44 

Nolin P and Werner W 2009 J. Am. Math. Soc. 22 797-819 

Platini T, Karevski D and Turban L 2007 J. Phys. A 40 1467-79 

Collura M, Karevski D and Turban L 2009 J. Stat. Mech. P08007 

Mandelbrot B B 1967 Science 156 636-8 

Delfino G, Viti J and Cardy J 2010 J. Phys. A 43 152001 

Nadarajah S 2005 J. Appl. Statist. 32 685-94 

Schramm O 2001 Electron. Commun. Prob. 6 115-20 

Simmons J J H, Kleban P, Dahlberg K and Ziff R M 2007 J. Stat. Mech. P06012 

Ikhlef Y and Ponsaing A 2012 iarXiv:1202.5476fc l 

Flores S M, Kleban P and Ziff R M 2012 larXiv: 1201 .64051 



